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Annotation 

Examined in the study is the probl<im of constructing the 
programmed motion of a multi-leg hounding vehicle in the flight 
phase. An algoritl.m is constructed for solving the boundary 
value problem of crnstructing the programmed motion for the 
flight phase. Also shovm is the possibility of the substantial 
simplification of the work of this algorithm for symmetrical 
motion of a vehicle. A method is proposed for constructing the 
translational motion of the legs, which provides no impact 
during lift-off and softness of placement of the legs on the 
supporting surface . The algorithms were developed by the 
method of mathematical simulation on a computer. The results 
of the calculations are given. An algorithm is realized for 
constructing the programmed motion utilizing tables for a 
sufficiently broad set of standard motion conditions. 
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Introduction 


The given study is a continuation of the investigations, 
begun in L33, on the creation of a mathematical model of motion 
and a motion control system for a bounding vehicle. Examined 
in study [3] was the supported phase of motion of the vehicle. 
Constructed in the present study is a mathematical model of 
motion of a vehicle in the flight phase, and the problem of 
construction of programmed motion in the unsupported phase 
is solved. 

A vehicle is examined which consists of a body and four 
or six two-section legs, which have three degrees of freedom 
each. The total weight of the legs comprises an appreciable 
portion of the weight of the body. Motion of the vehicle 
consists of the alternation of two phases: supported, when 
all of the vehicle's legs stand on a supporting surface and 
quasistatic stability ta!kes place, and unsupported, or the 
flight phase, in the course of which the center of mass moves 
along a ballistic trajectory. We request that the vehicle 
leave the supporting surface without impact, and, at the moment 
of touchdown, it should set down sufficiently softly, with a 
small absolute stop velocity (zero, if possible). The re- 
quirement for softness of placement of the legs, instead^ of 
impactless, ensures reliability of contact of the legs with 
the supporting surface at the moment of touchdown, with the 
presence of errors in execution and informational errors on 
the actual supporting area. 

In the problem of constructing the programmed motion of 
the vehicle in the flight phase, it is necessary to find those 
values of the linear and angular velocity of the yehicle 
which provide movement from a given initial position to a 
given final position, with the adopted method of shifting 
of the legs. 

Obtained in § 1 are the primary integrals of the equations 
of motion of the vehicle (law of motion of the vehicle's center 
of mass and law of conservation of angular momentum relative to 
the center of mass) , and a mathematical model of the spatial 
motion of the vehicle in the flight phase is constructed. 
Formulated in § 2 is the boundary value problem for construc- 
ting the programmed motion of a bounding vehicle in the flight 
phase, and an economical algorithm is obtained for its solution 
A method is proposed for constructing the translational motion 
of the legs, which ensures non- impact of lift-off and the re- 
quired softness of placement of the legs on the supporting sur- 
face. Also examined is the possibility of substantial simpli- 
fication of the working of the algorithm for constructing the 
programmed motion in the case of symmetrical motion of the ve- 
hicle. The results of the calculations are given in § 3 • 

For standard motion conditions, v/hich encompass a broad 
class of the most frequently encountered motion conditions of 


the vehicle, at is proposed that one recall the dependence of 
the initial angular velocity of the body on the parameters of 
the impending flight phase. The volume of stored information 
is relatively small. The value of the initial linear velocity 
of the body is calculated. Such an approach sharply shortens 
the time for calculations, and makes it possible to solve the 
problem of constructing the programmed motion for the flight 
phase during motion of the vehicle. 

The algorithm for constructing the programmed motion and 
the mathematical model of spatial motion of the bounding ve- 
hicle in the flight phase were realized in FORTRAN language on 
the BESI-6 computer. 

The constructed programmed motion can serve in the capacity of 
supported motion for construction of the algorithm of stabili- 
zation of the angular motion of the body during flight. 

The author expresses his deep appreciation to D. Ye. 
Okhotsimskiy for posing of the problem and attention to the 
study. 


§ 1. Mathematical Model of Motion of Bounding Vehicle in the /? 

Flight Phase 

1.1. Kinematics of Vehicle and Designations 

We will examine a vehicle which consists of a body and 
N(N=^,6) two-section legs. The motion of the vehicle occurs 
in the absolute system of coordinates The coordinate 

system 0^ is closely associated with the body of the vehicle, 
and its axes are directed along the main axes of inertia of the 
body. The position of the body of the vehicle we will deter- 
mine by the coordinates of the center of mass of the body 

and the angles 4',0 ,y (fig- 1) » where 4' is the angle be- 
tween the axis 0^^ and the projection of Oy on the plane Oj^n 
(yaw angle) , 0 is the angle between the axis Oy and its pro- 
jection on the plane Oj^n (pitch angle), and y is the angle 
of rotation around the axis Oy (bank angle) . The plane Ojfiri 
is horizontal. 

Let each leg consist of two sections (thigh and shank), 
and have two degrees of freedom at the point of support of the 
thigh and body, and one degree of freedom at the knee. The 
angle between the axis 0^ and the projection of the leg on 
the plane O^^y will be designated by a , the angle between the 
negative direction of the axis 0^, and the thigh by p , and the 
angle between the thigh and the shank by q (fig. 2). The plane 
of the leg, formed by the thigh and the shank, is perpendicular 
to the plane O^^y . 

The contact of the leg with the supporting surface has a 
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point nature. We will call the end of the leg the foot. The 
center of mass of the thigh lies on the axis which connects the 
joints at the point of support leg and body and at the knee t 
while the center of mass of the shank lies on the axis which 
connects the knee joint to the foot. 

We will introduce the following designations! 


M - mass of the body, 


'^x'Ty - main moments of inertia of the body, 


mj^ - mass of the j section of the j. leg (here, and subsequently 
j=l corresponds to the shank, and j=2 to the thigh) , 





tensor of inertia of 
the j. section of the 
i legi 


l®^S'th of the j section of the i leg, 

distance from the point of support of the j. section of 
the i leg to its center of mass. 

n »ynl> 2nl) " radius-vector of the point of support of 

the j. leg and body. 


fi'^^=(x^^- ,y^ 4 ,2^4) - radius-vector of the center of mass of the 
^ i section of the i leg. 


^ci“^^ci’^cl’ ^ci^ " radius-vector of the foot of the i leg, 

.4 - matrix of transition from absolute system of coordinates 
0i § n? to system of coordinates associated with the 

boay of the vehicle. 


p=(xjj,,yjvj, Zj^) - vector directed from the center of mass of the 
vehicle towards the center of mass of the body, 

R=(^,n,^) - absolute coordinates of the center of mass of the 
body 


I 1 - absolute coordinates of the center of mass of 

a. ^3. ‘a ’^ci .- . . - 

the vehicle, 


to - angular velocity of the body of the vehicle. 

We will write out some kinematic relationships: 



ifc; = , 


(I.I) 




xo =r^i*tC^i»p/^i Ai»^i)^i, 




^ii =ik- 

^,; = At; P <^A» 




/lO 


where 


-tt. 


^£. Mr, &.• 

f=- i-> “• J -i- 
M->-x~d^ **>/{ 

4^i i*i 


i- 


~A<ri'('t^& 




41ft 


( 1 . 2 ) 


1(1.3) 


1.2. Primary Integrals of Equations of Motion in the Flight Phase 

We will examine the motion of the vehicle with weighahle 
legs. During the flight phase, a single external force acts 
on the vehicle — the force of gravity, directed along the axis 
Oi^ and applied at the center of mass of the vehicle. 

From the theorem on motion of the center of mass of a 
mechanical system, it follows that the center of mass of the 
vehicle moves along a ballistic trajectory [2] 

where Ra“^“ radius-vector of the center of mass of the 
vehicle , 
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V„=R-5)Xp-p - velocity of the center of mass of the 
vehicle , 

w - angular velocity of the body. 

According to the theorem on the change in the vector of /ll 
the kinetic moment of a mechanical system, relative to the 
center of mass, we obtain [2] 


Qa. = , ^ 1 * 5 ) 


where G - vector of the kinetic moment of the vehicle in the 
absolute system of coordinates OiSnC. 

According to Koenig's theorem L23, the value of the vector 
of kinetic moment of a random solid body is equal to 


where t m ~ tensor of inertia of the body relative to its 
center of mass, 

- mass of the body, 

- radius-vector of its center of mass, 

- velocity of the center of mass, 

CO angular velocity of the body. 

We will make use of formula (1.6) for the calculation of 
the values of the kinetic moment of the body and the sections 
of the legs relative to the center of mass of the vehicle. 

As a result, we obtain the following expression for the 
vector of the kinetic moment of the body: 


+ M^(s, §) - Mf(c5,§) ( 1 . 7 ) 


where x - tensor of inertia of the body. 

We will associate a^ system of coordinates, designating its 
unit vectors » ejj_, with the j section of the i leg. 

The vector ej^'^i's directed along this section of the leg, and 
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,.,2 -3 , _1 -2 

eji is perpendicular to the plane of the leg, x ejiJ . 

The origin of the system of coordinatt^s ^ will he placed at the 
center of mass of the j. section of the 1 leg. In _ this case, 
the matrix of transition from the system of coordinates _0xyz 
to the system of coordinates associated with the j section of 
the i leg has the form [1] 


Jj. r I «>&<»{. 


coi<P-;Aind,' 



( 1 . 8 ) 


where 




y=y. 


The angular velocity of the j. section of the i leg, in 
the projection o'n the axes associated with this section, is 
equal to [1] 




Ca3 / 



( 1 . 9 ) 


The velocity of the center of mass of this section we will 
calculate according to Euler's formula [2] 


Vj; = f ^ i- 5v). (1. 10) 


Utilizing (1.6), (1.9), and (I.IO), we obtain the final ex- 


pression for the vector of the kinetic moment of the j, section 
of the i relative to the center of mass of the vehicle 


1 




(I. II) 


The vector of the kinetic moment of the vehicle, relative to 
the center of mass, is equal to the sum of the kinetic moments 
of the body and the sections of the legs 


Prom (l.l)-(1.3) and (1 .6) -( 1 .12) , it follows that 


- 1 ^ * 

/ p- 

where Pi=(a i,p i,qi) 

T 

f 


/hr 


4 ,\ 

/= / hjf 

4 » 


\hi 

h ' 

hi j 


(I.I3) 



Ijrfi. 








A 3 


In that case when the legs of the vehicle are immobile, relative 
to the body,_the kinetic moment of the vehicle is equal to the 
first term Iw in formula ( 1 . 13 ), and the second term of ^ 

is associated with motion of the legs relative to the body. The 
matrix I is the tensor of inertia of the vehicle with a given 
position of the legs in a relative system of coordinates. The 
elements of the matrices I,Ii (i=l,2,...N) depend solely on the 
angles in the leg joints ai, Pi, qi ( i=l , 2 , . . . N) . The values 
of the elements I,Ii are written out in explicit form in the 
appendix . 

Let control of shifting of the legs in the flight phase 
be accomplished according to position, i . e . . ai,Pi,qi (i = 
1,2,...N) are prescribed by functions of time. We will con- 
struct a mathematical model of motion of the vehicle, utilizing 
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the theorem on the motion of the center of mass and the 
theorem on the change in the vector of the kinetic moment 
relative to the center of mass. The equations of motion of 
the vehicle in projection onto absolute coordinates will be 
written in the foll^ '.ng manner [2] 


= a.w) 

( 1 . 15 ) 


In view of the nonsingular it;/ of the matrix I (I is the 
tensor of inertia of the vehicle), it follows from (1.5) i 

(1.13)1 and (1.15) that 



The right-hand portions of the equations ( 1 . 1 6) are functions 
of time and the angular coordinates of the body, since the A 4 

elements of the matrices I,!;}, depend solely on the angles 
in the leg joints ai, Pi*^ 3 ( i=l » 2 , . . .N) , which are known 
functions of time. 

We will add the kinematic equations of Euler [ 1 , 5 D to 

(1.16) 


( 1 . 17 ) 


By integrating equations ( 1 .16) -( 1 . 17 ) on a computer by 
the Runge-Kutt method, we obtain motion of the body of the 
vehicle near the center of mass, i . e . . the dependence of the 
angles b,0,Y on time. 

We would note that equations ( 1 . 1 ?) deteriorate with 
6=Tr/2. Geometric deterioration takes place. In that case 
when this may take place, it is sufficient to change the 
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direction of the axes in an appropriate manner. 

Motion of the center of mass of the vehicle will he de- 
termined from equations determined for a parabolic 

trajectory. 

The motion of the center of mass of the body will be 
found, utilizing the relationships* 




( 1 . 18 ) 


In the given paragraph, we first obtained the integrals 
of the equations of motion of a bounding vehicle in the flight 
phase. Based on the primary integrals, we constructed a 
mathematical model of motion of the vehicle in the flight 
phase . 


S 2. Construction of Prop:rammed Motion of the Vehicle for /l5 

the Flight Phase 

2.1. Posing of Problem 

The system of construction of the programmed motion of 
the vehicle in the flight phase determines the velocities of 
the body of the vehicle at the moment of lift-off from the 
supporting surface, necessary for realization of the impending 
flight phase, and creates the nominal translational motion of 
the legs. The initial and final position of the vehicle, 
determined by the algorithm of formation of the attitude of 
the bounding vehicle at the moments of lift-off and touchdown 
[_ 3 '], are the input parameters of the algorithm of construction 
of the programmed motion. In the flight phase, the vehicle 
should move so as to satisfy the boundary conditions: 


ja )= Y , A/). 


Where‘S =(C , rif V'jQ »y) 5 moments 
of beginning and ending of the flight phase. -^Positioning and 
lift-off from the supporting surface of all of the legs of the 
vehicle are accomplished simultaneously. 

In order to determine the initial velocity of the body of 
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the vehicl.e at the moment of lift-off from the supporting 
surface f(L)r i necessary for realization of the impending 
flight phase, knd to construct the nominal translational 
motion of the legs, we will utilize the following calculation 
scheme: 


1. Utilizing the houndary conditions (2.1), we will ob- 
tain the value of the velocity of the center of mass at the 
moment of lift-off from the supporting surface and the duration 
of the flight phase from the law of motion of the center of 
mass (1.4) . 

2. The translational motion of the legs in the flight 

phase leads to a change in the angular coordinates of the /l 6 

body, in view of the law of conservation of the vector of 

the kinetic moment (1.5). The problem of determining the 

initial angular velocity of the body amounts to a boundary 

value problem, which is solved by the iteration method. In 

each step of this iteration process, it is necessary to know 

the corresponding translational motion of the legs. 

The velocity of the center of mass of the body of the 
vehicle is determined according to Euler's formula [2], 
utilizing the known values of the velocity of the center of 
mass of the vehicle and the angular velocity of the body. 

3. Shifting of the legs will be accomplished in a rela- 
tive system of coordinates Oxy^, thus ensuring non- impact 

of lift-off and sufficiently soft placement of the legs on 
the supporting surface with a small absolute velocity, which 
comprises a given part of the velocity of the center of mass 
of the vehicle at this moment. At the moment of touchdown, 
the requirements for softness of placement of the legs, in- 
stead of non- impact, ensures reliability of contact of the 
legs of the vehicle with the supporting surface in real 
motion, with the presence of errors in execution and infor- 
mational errors on the real support area. 

The boundary values of the velocities at the leg joints 
depend on the initial velocity of the body of the vehicle. 
Consequently, construction of the translational motion of 
the legs must be accomplished in each step of the iteration 
process of the solution of the boundary value problem by 
determining the initial angular velocity of the body. 

We will switch to the description of the operation of 
the basic blocks of the algorithm of construction of the 
programmed motion. 


2.2. Determination of Duration of Flight Phase and Parameters /± 7 

of Trajectory of Motion of Vehicle's Center of Mass 
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From the houndary conditions (2.1), which determine the 
initial and final position of the body and feet of the vehicle 
in the absolute system of coordinates Oi^nC» we will find the 
matrix .4 , the relative coordinates of the feet and angles in 
the leg joints, and the value of the vector p", directed from 
the center of mass of the vehicle towards the center of mass 
of the body, at the moments of lift-off and touchdown. We will 
then obtain the values of the coordinates of the center of mass 
of the vehicle at the moments of lift-off and touchdown 




Let the tangent of the angle of inclination of the 
initial velocity of the center of mass of the vehicle to 
the horizontal plane Oj^n, equal to X, be given. 

Prom the law of motion of the center of mass of the 
vehicle (1.4) and the boundary conditions (2.2), we will 
find the duration of the unsupported phase of motion and the 
initial velocity of the center of mass of the vehicle 


v.;-a‘-0/T , 


( 2 . 3 ) 


Equations (1.4) and (2.3) fully determine the motion of 
the center of mass of the vehicle in the flight phase. 

2.3. Shifting of the Legs 

According to the absolute coordinates of the feet and body /18 
of the vehicle, known, from the boundary conditions (2.1), we 
will calculate the coordinates of the feet in the relative 
system of coordinates and, transforming the third group 

of equations (1.1), we will obtain the values of the angles at 
the leg joints at the moments of touchdown and lift-off from 
the supporting surface, according to the formulas: 
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Formulas (2.4) deteriorate when the projection of the foot 
on the plane Oxy coincides with the projection of the point 
of support of the leg on this plane, i .e . , ^jt’sin ®x+ j? j j. 
AjfsinPi=0. With the formation of the inixial and final 
attitude of the vehicle, we require that, at these moments 
in time, the projections of the foot and the point of support 
of the leg on the plane Oxy do not coincide. 

The values of the velocities at the leg joints at the 
moments of lift-off and touchdown are determined from the 
conditions of non-impact of lift-off of the legs from the 
supporting surface and the conditions of sufficiently soft 
placement of the legs with a small absolute velocity (zero, 
if possible) , which comprises a given part of the velocity of 
the center of mass of the body at the moment of touchdown 




•r 

f 


( 2 . 5 ) 


where Cg,C^,C;-«l are the given constants. 

The velocity of the center of mass of the vehicle is known 
from (1.4) and (2.3); thus, the absolute velocity of the feet 
at the moments of lift-off and touchdown are known. We will 
write out the expressions for the absolute velocity of the 
feet, according to Euler's formula [2 3 

a ^ V- J , (2.6) 


We will substitute in (2.6) the values of p, n cl ( i=l » 2 , . . .N) , 
obtained by differentiation of (1.1) -(1.2); then, (2.6) form 
a system of linear equations relative toai,0i,qx ( i=l , 2 , . . .N) , 
which we will write in the form 


/ s =■ _ ( 2 . 7 ) 

^ P~‘ S ^ k-fZ, ^ , 


where p=( , P q. q 2 > • • ) “3N is the measured vector of the 

angles at the leg joints, L is the vector of dimensionality 


3N X 3N, the elements of which depend solely on p; S=(7a-Vcl» 

. . .V^-Vcn) { P, is _the_matrix of dimensionality 3 x 3N, so that 
(5)X ( p’^n 0 ^ * ' ' ' ^ c^i^ ' 

As a result, we obtain the fact that, having determined 
the angles at the leg joints and constructed the trajectory 
of motion of the center of mass of the vehicle, one can deter- 
mine the values of the elements of the matrix L, f2 and the 
vector S at the moments of lift-off and touchdown. If, in 
addition, the initial and final value of the angular velocity 
of the body are known, then,j|_by solving the systems of linear 
equations (2.7) relative to p, we will find the boundary 
values of the velocities at the leg joints at the moments of 
lift-off and touchdown. The system of equations (2.?) 
characterizes the dependence of the velocities at the leg 
joints on the velocity of the center of mass of the vehicle 
and the angular velocity of the body in that case when the 
coordinates of the body of the vehicle, the angles at the leg 
joints, and the absolute velocities of the feet are known. The 
matrix L in (2.?) deteriorates when the velocity of the foot, 
as if it were one of the legs, does not determine the velocity 
of the joints of this leg unequivocally in the relative system 
of coordinates. This takes place when one of "the following 
conditions is fulfilled: a) the projection of the foot of the 
i, leg on the plane coincides with the projection of the point 
of support of the leg and body, or b) the leg folded up or 
completely straightened at the knee (qi=0 or qi=ir)._ During 
formation of the initial and final programmed positions of the 
vehicle, we require that these conditions not be fulfilled. 
Then, the matrix L is n*'t deteriorated. 

With each return to the block of construction of the 
translational motion of the legs, the value of the initial 
angular velocity of the body u ( tp ) = u)° is known. The value 
of the angular velocity of the body at the moment of touch- 
down will be determined utilizing the law of conservation of 
the vector of the kinetic moment. Prom (1.5) and (1.13), it 
follows that 


^Jl f) = dj" ( 2 . 8 ) 

T 

where Il = (Ij_ ,l 2 » • • . % ) • 

By substituting the values of p(%) , obtained by solving the 
system of linear equations (2.7) and (2.8), we obtairua 
system of three linear equations relative to w(tT)= 

= C-rV, L'S, ^ 2 . 9 ) 
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by solving which we will find w' • 

As a result, we obtain the values of the coordinates and 
the velocities at the_ leg joints at the moments of lift-off 
and touchdown. Let p=( a^, Pi, . . .q[^- ) • Shifting of the legs in 
the flight phase will be accomplished so that the obtained 
boundary conditions will be fulfilled: 


paj= F* , p(L)--p^ 


(2.10) 


and the phasic restrictions will not be disrupted: 




(2.II) 


where given limits of change- in the angles 

at the leg joints in the programmed motion. 

Then, for Pi ( i=l , 2 , . . .N) , we have 


, 'PUi- i. ) » A" 

p->o , pUirU »pt-pr. 

* 


Z21 


Consequently, the shifting of the legs should be accomplished 
so that Pi(t) fl'>"st slows from the initial position on the 
phasic plane (PiiP^) to (Pj^^’^O) , remains for some time at 
the boundary (p^min,0), and then accelerates to a final position 
(0-^.Pi’*) • Motion at the knee is accomplished similarly, only 
qi ( c) (i=l,2...N) reaches the upper limit In other 

words, right after lift-off from the supporting surface, the 
leg straightens and increases its moment of inertia relative 
to the axis of the engine, which accomplishes a change in the 
angle Pi . For more uniform distribution of the forces developed 
in the leg joints, according to time, we require that the 
accelerations along the coordinates p. ,q^ , with braking of the 
leg after lift-off from the supporting surface, be linear 
functions, which diminish, according to the modulus, to a 
zero value. In a symmetrical manner, we v/ill accomplish 
motion along the coordinates pj_ ,qj_ right after_ the moment of 
touchdown. We will designate t = Pi, q i, PiCiin^ ( i=l , 2 

. . .N) . We will require that the boundary conditions (2.10) be 
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fulfilled and 




/* a J ^ » u/iTri 

(2.12) 



Uiho , UCTXA. 

(2.13) 


Then, it follows from (2.10), ( 2 . 12) -( 2 . 13) that 


h - to 
r j 

' gt 




a, 


(2.R) 


From (2.12)-(2.14) , we obtain the dependence on time of the 
coordinates Zj=Pi,qi (i=l,2...N) with shifting of the legs 
in the flight phase 


% Q^a-QY3ti -Ho -O , rA 

ft'’'’ . U.CTH (2. 15) 


At the same time, for the engine, which accomplishes 
turning of the plane of the leg (i=l,2...N), the moment 
of inertia of the leg changes slightly. Therefore, according 
to the coordinates , shifting of the leg will be accomplished 
in such a manner that will minimize the maximum magnitude of 
acceleration according to this coordinate, satisfying the boun- 
dary conditions (2.10) and the phasic limits (2.11). Con- 
struction of motion along the coordinate a:t amounts to the 
solution of the extremal problem 


( 2 . 16 ) 


4 ' 


A" <x(/ , Ji‘ fL)~^i , 


Problem (2.16) was examined in study L33» where its 
analytic solution was obtained. 

As a result, we obtained the dependence of the angles at 
the leg joints in the fliglit phase on time and the initial 
angular velocity of the body. 


2.4. Boundary Value Problem. Determination of Initial Angular /23 

and Linear Velocity of Body. 

We will formulate the boundary value problem for deter- 
mining the initial angular velocity of the body. By substi- 
tuting the dependence, obtained in paragraph 2.3, of the 
angles at the leg joints on time and the initial angular 
velocity of the body p=p(t,w% where P=( ai, » • • • ,q.N) , 

into the law of conservation of the vector of -fehe kinetic 
moment (1.5) i ( 1 . 13 ) » we have 


Kiroo)-. ^ 5/4 


(2. IV) 


Prom the kinematic equations of Euler (1.17) and equations 
(2.17), it follows 
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IfJe will write out the boundary conditions (2,1) for the angular 
coordinates of the bodys 


mo) = r , eao)= o'" , ^(^o) ! 


We will solve the boundary value problem ( 2 . 1?) -( 2 , 19) for 
the system of differential equations, which describe the 
motion of the vehicle near the center of mass in the flight 
phase, numerically on the computer. As a result, we obtain 
the value of the initial angular velocity of the body or® , 

By integrating system (2.18) in the time increment 
■fceCtC'tJ (where to is the beginning and ti is the end of 
the flight phase) by the Runge-Kutt method, we obtain the 
values of the angular coordina'Jies of the hOpdy at the moment 
of touchdown as functions of or ij/ , 6 “, y The boundary 
value problem (2.17)-(2.19) comes down to the solution of 
the system of algebraic equations 




The system of equations (2.10) was solved numerically on the 
computer by the modified method of the quickest ascent [4,63* 

We would note that translational motion of the legs is 
created in each step of the iteration process of solution of 
the system of equations (2.10), and, consequently, the boundary 
values of the velocities at the leg joints are determined. For 
this purpose, it is necessary to twice solve the system 3N of 
linear equations (2.7). The volume of calculations in each 
step of the iteration process is substantially reduced if, 
prior to solving the boundary value problems, one solves the 
corresponding matrix equations by the Gauss method and finds 

and L”^ at the moments of lift-off and touchdown. In 
this case, in each step of the iteration process of solving 
the boundary value problem, instead of solving the two systems 
3N of linear equations (2.7), it is sufficient to calculate the 
boundary values of the velocities at the leg joints, according 
to the formula 


C2.2I) 


p-fl-S)i-(L-^)^ It./,, - • 


that is, it is sufficient to multiply the matrix L"* by the 
value of the vector of the angular velocity of the body, and 
to add it to the vector L”^5. 

Knowing the angular velocity of the body at the moment /25 
of lift-off, as a result of solving the boundary value prob- 
lem (2.17)-(2.19) » and the initial velocity of the center of 
mass of the vehicle from (2.3), we will find the velocity of 
the center of mass of the body at the moment of lift-off, 
according to Euler's formula 12] 

( 2 . 22 ) 


2.'^. Simnlified Scheme of Construction of Programmed Motion 

Let the vehicle have an even number of legs N=2n, sus- 
pended symmetrically relative to the plane Oyz» and the right 
and left leg of each pair of symmetrically suspended legs have 
identical dynamic and kinematic characteristics. 

Let the positions of the vehicle at the moments of lift- 
off and touchdown satisfy the conditions: 


nio)-'f'Ch) , (2.23) 


and the right and left leg of each pair of legs be located 
symmetrically; then, in the flight phase, the right and left 
legs will move symmetrically, and the angles of yaw and bank 
will not change 


( 2 . 2 ^) 


In this case, the problem of constructing the programmed 
motion can be simplified substantially. Taking into account 
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the symmetry of motion, the law of conservation of the kinetic 
moment of the vehicle, relative to the center of mass, can be 
written in the form 


where 






' fh <■ A *1^- ijr. (2.j;5) 


C/cu s O " O 


, ^/= J=^>~ 


4^. = y 

^ , /p ,/^/. 


We will solve the boundary value problem according to the 
determination of the initial angular velocity of the body, 
according to the single variable 0 A For the determination 
of the initial and final velocity at the leg joints, it is 
sufficient to solve the two systems 3N/2 of linear equations 
instead of the two systems 3N of linear equations (2.7), that 
is, lesser dimensionality by twice than with the full scheme 
of calculations. All the rest of the calculations are also 
substantially simplified. With regard for the symmetry of 
motion of the bounding vehicle in the flight phase, planar 
motion of the vehicle is examined, instead of spatial motion, 
during construction of the programmed motion according to the 
simplified scheme. 

2.6. Algorithm for Constructing the Programmed Motion Utilizing 

the Solution of the Boundary Value Problem 

We will examine the working of the algorithm for con- 
structing the programmed motion of the vehicle in the un- 
supported phase of the bound as a whole. We will show the 
sequence of the calculations. The block diagram of the algo- 
rithm is given in figure 3; in this case, the number of blocks 


corresponds to the numbers of points in the description of 
working of the algorithm given below. The algorithm auto- 
matically creates a selection of the simplified or complete 
scheme of calculations, carried out in blocks 3-6. 

1. We will calculate the values of the matrices of tran- 
sition (1.3), the relative coordinates of the f^et and the 
angles at the leg joints (2.?), and the vector p, directed 
from the center of mass of the vehicle towards the center of 
mass of the body (1.1) -(1.2), according to the absolute coordi- 
nates of the f&et and the coordinates of the body of the vehicle, 
known from the boundary conditions (2.1), at the moments of 
lift-off from the supporting surface and touchdown. 

2. We will determine the absolute coordinates of the 
center of mass of the vehicle at the moments of lift-off and 
touchdown (2.2). We will calculate the duration of the 
flight phase and the initial velocity of the center of mass 
of the vehicle (2.3). As a result, the ballistic trajectory 
of motion of the center of mass of tl',. vehicle is fully deter- 
mined . 


3. We will calculate the elements of the matrix L, and 

the vector 7 , which are part of equations ( 2 . 1 ^), at the 
moments of lift-off and touchdown. By solving the _ corresponding 
(2.14) matrix equations , _by the Gauss method, _we will find the 
elements of the matrix and the vector L” 5 with t=to and 

t=t]_ (here and subsequently, tq is the moment of the beginning 
and ■tl is the moment of completion of the flight phase) . 

4. We will find the initial angular velocity of the body 
of the vehicle, and the boundary value problem (2.18)-(2.19) » 

which amounts to the solution of the system of equations ( 2 . 20 ). /29 

Then, having calculated the initial velocity of the center of 
mass of the body ( 2 . 22 ), the algorithm completes the work. 

Calculation of the left parts of equations (2.20), in the 
process of the solution of the boundary value problem, takes 
place in blocks 5»6. We would note that, with each treatment 
of these blocks, the initial angular velocity of the body of 
the vehicle is known. 

5 . Also determined is the angular velocity at the leg 
joints at the moment of lift-off (2.21) , and the vector of 
the kinetic moment of the vehicle (I.I3) is calculated. By 
solving the system of three linear equations (2.16) , we will 
determine the angular velocity of the body at the moment of 
touchdown. From relationship (2.21), written out with t=tj , we 
will find the velocity at the leg joints at the moment of touch- 
down. We will construct the translational motion of the legs. 

6. By integrating the system of differential equations of 
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Fig. 3 


Key: a. Input 

to. Determination of charac- 
teristics of initial and 
final attitude of vehicle 

c. Construction of trajec- 
tories of center of mass 
of vehicle 

d. Preliminary calculations 
for determining the veloci- 
ties at the leg joints at 
the moments of lift-off and 
touchdown 


e. Solution of the 
tooundary value protolem 

f. Construction of the 
trajectories of shifting 
of the legs 

g. Integration of the sys- 
tem of the kinetic 
moment 

h. Output 


the primary integral — the law of conservation of the vector 
of the kinetic moment (2.18) toy the Rungt,-Kutt method, with 
given initial conditions and the const^-ucted shifting motion 
of the legs, we v/ill find the corresponding given initial con- 
ditions of the value of the angular coordinates of the toody at 
the moment of touchdown. 

Note . In blocks 3-6, all of the references pertain to the 
formulas which correspond to the complete scheme of calculations. 
In the case of working of the algorithm according to the simpli- 
fied scheme (paragraph 2.5) » they are simplified, with regard 
for the symmetry of motion which takes place. 


3^ Construction of Programmed Motion of Vehicle Utilizing 

Tables. Standard Motion Conditions. Results of Calculations 



Calculations were carried out for a four-leg:ged^ vehicle , 
having the following dynamic and kinematic characteristics • 

The mass of the "body is 1,000 kg, and the main moments of 
intertia (in kgm ) are: Xx =530i''^y = 90 , xs =605 • A12. of the 
legs are identical, and consist of two sections (thigh and 
shank). The length of the thigh is 1.2 m, and the length of 
the shank is 1.^ m. The center of mass of each section is 
located in the middle. The mass of the thigh is 30 kg, the 
mass of the shank is 20 kg, the main moments of inertia of 
the thigh (in kgm ) are: Xxi =3*2, Xy 2 =0.05, Xs 2 =3*2, and 
those of the shank are: xa:]_= 3 - 6 , tyx= 0 . 05 , X 2 x“ 3 -^'* The 
coordinates of the points of support of the legs and the 
body have the values: Ki =0.5 ( 1 = 1 , 2 ), x =0.5 (i=3»^)» 2/1= 

1.25 (i=l,3), (1=2,4), ^ 1=0 ( 1 = 1 , 2 , 3, 4). 

The system of units of measurement is the International 
System of Units. The time is measured in seconds, and the 
angles in radians. Acceleration of the force of gravity 
^= 9*8 m/sec. 

The initial and final programmed positions of the vehicle 
are determined by the algorithm of formation of the attitude 
of the vehicle at the moments of lift-off from the supporting 
surface and touchdown, the working of which is examined in 
detail in [3]. The index contour in the projection on the _ 
horizontal plane Ox^n is a rectangle, the length of v/hich is 
equal to the distance between the points of support of the 
front and rear legs, and the width is equal to twice the 
lateral extension of the legs from the longitudinal axis of 
"he index contour, which we will designate using 7^* The 
input parameters of the algorithm of formation of the attitude, 
the selection of the values of which is accomplished either by 
the operator of the bounding vehicle or by higher levels of 
the motion control system, are: D - the distance of the bound 
(distance between the centers of the index contours in the pro- 
jection on the horizontal plane Ox?n) , 'jJf - the direction (azi- 
muth) of motion in the flight phase, the orientation 

(azimuth) of the initial index contour, ip'ftji - the orientation 
(azimuth) of the final index contour, X - the tangent of the 
angle of Inclination to the horizontal of the initial velocity 
of the center of mass of the vehicle, and also information on 
the supporting surface. At the moments of lift-off from the 
supporting surface and touchdown, the distance along the ver- 
tical from the supporting area to the center of mass of the 
body is equal to the given value Cnom , the angles of pitch and 
bank are equal to zero , 0 '’= y®= y '''=0 > the angle of yaw^ is 

equal to the azimuth of the index contour, ='PuA» -’I'ttA . 

The coordinates ^ and n of the center of mass of the body at 
the initial and final moments are determined from the condition 
that, in the projection on the horizontal plane, these points 
lie on a line connecting the__ centers of the index contours at a 
distance ( ?ncm - Cmir) /ii^ax( 1 , \* ) from the center of the corres- 
ponding index contour. Here, ^min is the minimally permissible 


distance along the vertical from the supporting surface to 
the center of mass of the body, and is the tangent of the 
angle of inclination of the velocity of the center of mass of 
the vehicle to the horizontal at the corresponding moment in 
time . 


For the vehicle being examined, we will set ?mln“®*5 in* 

In study [3]; it was shown that the best values, from the point 
of view of minimization of the maximum magnitudes of the 
forces developing in the leg joints in the supported phase 
of motion, are the values ?nom =i*7 m, Cniln"^ Subsequently, 
we will consider that Cnom take oh these values. 

The algorithm of construction of the programmed motion 
of the bounding vehicle was realized in FORTRAN language on 
the BESI-6 computer. Translation was carried out by the 
FOREX optimizing translator, developed in the Institute of 
Applied Mathematics of the Academy of Sciences of the Soviet 
Union. Construction of the programmed motion of a four-legged 
vehicle in the flight phase, according to a simplified scheme 
of calculations for symmetrical motion of the vehicle, takes 
0.58 seconds, and for the complete scheme of calculations, 2.9 
seconds are required. On the Earth's surface, the duration of 
the unsupported phase of the bound, a distance of 10 m along /32 

the horizontal plane, is roughly 1 second. Therefore, it 
is necessary either to calculate the programmed motion for the 
bounding vehicle in advance for the series of bounds, or to 
have the solution of the problem of construction of the pro- 
grammed motion for the plight phase in the memory of the com- 
puter, which controls the motion of the vehicle. It is 
advisable to utilize both of these methods for constructing 
the programmed motion. For a set of standard motion con- 
ditions, encompassing a sufficiently broad class of the most 
frequently encountered motion conditions, it makes sense to 
retain the solution of the problem of constructing the 
programmed motion for the flight phase in the m^emory of the 
computer. For realization of a series of non-standard bounds, 
it is necessary to stop the vehicle and make preliminary cal- 
culations of the programmed motion. 

The basic volume of calculations during construction of 
the programmed motion of the vehicle for the unsupported phase 
of the bound is associated with the solution of the boundary 
value problem on the determination of the initial angular 
velocity of the body. For purposes of economy of memory of 
the controlling computer, necessary for remembering the sol- 
ution of the problem of constructing the programimed motion in 
the flight phase for standard motion conditions, it is advisable 
to store only the magnitude of the initial angular velocity of 
the body. The initial velocity of the center of mass of the 
body is calculated according to the finite formulas (2.3), 

( 2 . 22 ). In this case, with a comparatively small volume of 


26 



stored information, the volume of calculations for construction 
of the programmed motion for the unsupported phase of the bound 
is reduced substantially. Construction of the programmed 
motion, utilizing the tables, for standard motion conditions 
makes it possible to reduce the time of solution of this prob- 
lem up to 0 .058 seconds ( or roughly by $0 times, as compares 
with the solution of the problem of construction the programmed 
motion utilizing the boundary value problem) for the complete 
scheme of calculations, and for the simplified scheme of cal- /33 
culations — up to O .037 seconds (roughly by I 5 times) . 

Analogous logic is utilized in that case when the vehicle 
calculates the programmed motion for the series of bounds in 
advance. Only the initial angular velocity of the body is 
stored. For realization of motion of the vehicle, the value 
of the initial velocity of the center of mass of the body is 
calculated prior to each bound, with regard for corrections 
which may occur in the process of m.otion (for example, for 
the range of the bound) . 

We will call turning in the flight phase correct if the 
conditions 'kf fulfilled for the azimuth angle, 

that is, the direction of motion of the vehicle in the flight 
phase is the bisector of the angle formed by the longitudinal 
axes of the index contours at the moments of lift-off from the 
supporting surface and touchdown. The magnitude 
will be called the angle of turning of the vehicle in the flight 
phase . 

We will turn to the description of the standard motion 
conditions of the vehicle. We require that each of the 
supporting areas differ little from the section of the hori- 
zontal plane. The locale, as a whole, may have a complex 
form. Under standard motion conditions, the vehicle executes 
only correct turns within the limits of 1 0 .6 radians. 

The range of the bound lies within the range of from 5 to I 5 m, 
or the distance between the front legs of the vehicle at the 
moment of lift-cff from the supporting surface and the rear 
legs at the moment of touchdown in the projection on the hori- 
zontal plane 0 i 5 h is from one to five times the body of the 
vehicle. The tangent of the angle of inclination of the 
initial velocity'- of the center of mass of the vehicle to the 
horizontal either takes on one of the standard values 


A =/o.6, 0.7, 0.8, 0.9,1. 0,1. 25. 1. 5, 1. 75,2.0,2.5], (3.1) 


or is calculated in such a manner that the modulus of velocity 
of the center of mass of the vehicle at the moment of lift-off /34 

from the supporting surface for the given bound is minimal. 
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The value of X <1 (low-angle trajectory) should be utilized in 
the "run" condition along a relatively smooth surface [ 3 ]* The 
value of \>1 (steep trajectory) should be utilized for bounds 
on a slippery surface, for bounding across an obstacle of the 
barrier type, and for overcoming high projections. 

The solution of the problem of constructing the programmed 
motion of the bounding vehicle in the flight phase, for stan- 
dard motion conditions, depends on the following magnitudes: 

D - range of the bound, 

- angle of turning of the vehicle in the flight 
phase , 


- difference in altitudes of thi initial and 
final supporting areas or the height of 
the projection being overcome, 

X - tangent of the angle of inclination of the initial 

velocity of the center of mass of the vehicle to the 
horizontal . 

According to the results of the calculations, we will con- 
struct a table of the dependence of the initial angular velocity 
of the body on the magnitudes of D,ACttA » 

will store the values of the initial angular velocity for the 
junction points of the table, and for the actual values of D, 

ASttAi and AiI^ttAi the magnitude of the initial angular velocity 
will be determined by the method of linear interpolation. The 
junction points are selected so as to ensure accuracy of inter- 
polation of 0.05 radians/second. 

The complete table of the dependence of the initial angular 
velocity of the body on D,A!;ttA ,Ai|>tva , and X for the indicated 
motion conditions has a sufficiently large volume (about 15,000 
elements) . For purposes of decreasing the volume of stored 
information, we will narrow down the class of standard motion 
conditions. We will exclude the seldom utilized possibilities, 
placing additional restiictions of the values of the para- 
meters X and A 4 ;^y^, as a function of the height of the pro- /35 

jections being surmounted. 

Let Dq be the distance between the points of support of 
the front and rear legs; then, Df=D-Do is approximately equal 
to the distance by which the center of mass of the vehicle 
is displaced in the flight phase. 

According to the height of the projections being sur- 
mounted in the given bound, we will divide them into three 
classes: small projection | A^ttA l'^ ^ • ^5 D_, moderate projection 
0.05 D -^ ! A^Try^ [^0.25 D , high pro jection^O . 25 I A^irA *5 
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We will examine the standard motion conditions which correspond 
to each class of height of the projection being surmounted. 

Small projection |A?irA!^0‘05 Df* The value of the initial 
angular velocity of the body, in this case, practically does 
not depend on 

Given in figure 4 is the dependence o,f the initial angular 
velocity, according to the angle of pitch 6® , on the distance 
of the bound D for some values of the tangent of the angle of 
inclination of the initial velocity of the center of mass of 
the vehicle to the horizontal X for bounds without turning 


In order to ensure the required accuracy of interpolation, 
it is sufficient, for each X, to store the value of the initial 
angular velocity for five junction values of D which are 
multiples of the length of the body 


J ) 1 /'5.0, 7. 5, 10,0,12.5,15.0/. (3.2) 


The tangent of the angle of inclination of the initial velocity 
of the center of mass of the vehicle to the horizontal takes 
on one of the standard values (3»1)> 

Given in figure 5 is the dependence of the initial angular 
velocity of the body on the angle of turning of the vehicle 
in azimuth in the flight phase with a bound distance D= 

10 m for some values of X. The values of the initial angular 
velocity of the body, obtained through calculations, approxi- 
mate the linear functions well. Analogous results of calcu- 
lations for other values of D and X show that it is sufficient 
to have three junction values of ^'t'lrA 


=/o.o,o.3,o,c/ 


(3.3) 


With a change in the sign of the angle of turning in 
azimuth in the flight phase Aij^irA » in view of the symmetry of 
the vehicle, the initial velocity, according to the angle of 
pitch, is not changed, while the velocities according to the 
angles of bank and yaw change signs. This makes it possible 
to determine the initial angular velocities of the body for 
negative angles of turning Ai}jt 7 A< 0 » having in the memory the 
values of the Initial angular velocity of the body which 
correspond only to positive angles of turning Atjj.j^yy>0. 





The stored table of values of the initial angular velocity 
consists of 350 elements. It is taken into account here that, 
with A(fi 7 TA= 0 , it is necessary to store only the values of the 
initial velocity according to the angle of bank 6® (initial 
velocities according to the angles of pitch and yaw, in this 
case, are equal to zero, ® =y®=0) . 

High projection 0 .25Df $ | ACttA .5Df • We require that, 
in surmounting a high projection, the vehicle does not execute 
a turn in the flight phase, according to the angle of yaw, i.e. , 
At|^ 7 rA= 0 . The tangent of the angle of inclination of the initial 
velocity of the center of mass of the vehicle to the horizontal 
will be selected so as to minimize the value of the modulus 
of velocity of the center of mass of the vehicle at the moment 
of lift-off from the supporting surface. Taking into account 
the fact that the center of mass of the vehicle in the flight 
phase moves along a parabolic trajectory (1.4), we will obtain 
the optimal value 



where is the distance traversed by the center of mass of 
the vehicle in the flight phase in the projection on the hori- 
zontal plane OjCn. In this case, the modulus of velocity of 
the center of mass of the vehicle at the moment of touchdown 
will also have a minimum value. Utilized in the calculations 
was an approximated formula, obtained by the substitution of 
Da for its approximated value Dp in (3*4). 

The dependence of the initial angular velocity of the 
body of the vehicle on the range of the bound for some values 
of the relative height of the surmounted projection 
is given in figure 6. 

It follows from the calculations that it is sufficient to 
have five junction values for the range of the bound (3*2) and 
twelve junction values of 


The table of corresponding values of the Initial angular 
velocity of the body contains 60 elements. 

Moderate projection 0.05Dp « | 0. 25Dp. In this case, 

we will examine two standard motion conditions: 

1. The vehicle may execute a turn in the flight phase 




according to the angle of yaw. The value of the tangent of 
the angle of inclination of the initial velocity of the center 
of mass of the vehicle we will calculate as well, just as for 
surmounting a high projection (3*^) • 

Given in figure 7 is the dependence of the initial 
angular velocity of the tody of the vehicle on the angle of 
turning in the flight phase with a range of the bound 

D=10 m, for some values of the relative height of the sur- 
mounted projection 6=AC^y^/D|>. 

It follows from the results of the calculations that it 
is sufficient to have five junction values of the range of 
the bound D (3*2), three junction values of the angle of 
turning AijiirA (3*3) » and six junction values of the difference 
in heights of the supporting areas ACttA 

( 3 . 6 ) 


The stored table, in this case, contains- 210 elements. /42 

2. The vehicle surmounts an average projection without 
turning according to the angle of yaw Aip7rA=0* The tangent 
of the angle of inclination of the initial velocity of the 
center of mass of the vehicle to the horizontal may take on 
one of ten standard values (3*1) • We require that, at the 
moment of touchdown, the tangent of the angle of inclination of 
the velocity of the center of mass of the vehicle to the hori- 
zontal is at least 0.5» Then, \'e obtain the stricter restric- 
tion for the maximum value ofA^TrA with 0 


where 









If a signal with the required values of X, D, and ^^ttA enters 
the input of the algorithm for constructing the pro|:rammed motion, 
so that ACtj/ ( 4 ,D^) , then the least value of X" is selected 



from (3*1) I with which • This value of the 

tangent of the angle of inclination of^the initial velocity of 
the center of mass of the vehicle to the horizontal is con- 
sidered programmed X=X*-. 

For the required accuracy of interpolation, it is sufficient 
to have five junction values of D (3*2) and six junction values 




os.df ,-e.Z 


( 3 . 8 ) 


For the given motion conditions, it is necessary to 
store a table of initial values of the angular velocity of 
the body consisting of 300 elements. 

As a result, we obtain the fact that for storing of the 
initial angular velocity of the body of the vehicle for all of 74-3 
the examined standard motion conditions, it is necessary to 
store a table of 920 elements. 

For standard motion conditions, the valur-*! of the initial 
angular velocity of the body is determined from tho table, 
stored in the memory of the control computer, and the initial 
velocity of the center of mass of the body is calculated 
according to finite formulas (2.3) » (2.22). Such an approach 
maJces it possible to solve the problem of constructing the 
programmed motion of the vehicle for the flight phase in the 
process of motion of a bounding vehicle. With a comparatively 
small volume of stored information (920 numbers) in the BESI-6 
computer, construction of the programmed motion takes O .058 
seconds for the complete scheme of calculations, and 0.037 
seconds for the simplified scheme of calculations. 

Similar results were obtained for a six-legged vehicle. 

In this case, construction of the programmed motion, utilizing 
the solution of the boundary value problem for determining the 
Initial angular velocity of the body, takes ^.6 seconds for 
the complete scheme of calculations, and 0.86 seconds for the 
simplified scheme of calculations. For standard motion con- 
ditions, construction of the programmed motion for the flight 
phase, utilizing the tables, takes 0.097 seconds for the com- 
plete scheme of calculations, and 0.054 seconds for the simpli- 
fied scheme of calculations. 

Working out of the algorithms of the system of con- 
struction of the programmed motion of the vehicle in the un- 
supported phase of the bound, by the method of mathematical 
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simulation on the computer, showed the effectiveness of the 
proposed solutions. 

We will study the question of the dependence of the 
solution of the problem of constructing the programmed motion 
of the vehicle for the unsupported phase o” the bound on the 
relationship of the mass of the body and the total mass^o:^ the 
legs. We will designate this relationship through \x=W^^ 

We v/ill establish the mass of the legs, and will change the 
mass of the body. We will examine bounds without turning of 
a four-legged vehicle on a horizontal plane at a distance D. 

The angle of inclination of the initial velocity of the 
center of mass of the vehicle to the horizontal is equal to 
45 °. In view of the symmetry of motion, the angles of yaw 
and bank are identically equal to zero. Given in figure 8 is 
the dependence of the initial angular velocity of the body, 
according to the angle of pitch 0 ®, on (i for some values of 
the bound range D. With the mass of the body is con- 

siderably greater than the total mass of the legs, and, conse- 
quently, the translational motion of^ the legs practically has 
no effect on the motion of the bodye®-0. With p-O, the initial 
angular velocity of the body appro§.ches The energy of 
rotational motion of the body |Tjj(eo )2 is a monotonously 
diminishing function of |i. That is, diminishes more rapidly 
than l/VfT . On the whole, for the examined cases of ji^0.5i 
the energy of rotational motion of the body comprises no more 
than 1% of the kinetic energy of motion of the center of mass 
of the vehicle at the moment of lift-off from the supporting 
surface . 

The total energy outlays in the supported phase of the 
bound from place are associated with the transmission of 
kinetic energy to the center of mass of the vehicle and the 
kinetic energy of rotational motion to the body and leg sec- 
tions by the moment of lift-off from the supporting surface, 
as well as with a change in the potential energy of the vehicle. 
We will designate the energy outlays in the supported phase of 
the bound from place using and the specific energy outlays 

per unit of mass of the vehicle using We will 

study the question of the dependence of the specific energy 
outlays ^ on the range of the bound without turning on the 
horizontal plane D and the relationship of the mass of the 
body to the total mass of the legs M' • The angle of inclination 
of the initial velocity of the center of mass of the vehicle 
to the horizontal is equal to 45°. For the examined cases of 
li^O.5, the specific energy outlays practically do not depend 
on (I. The dependence of the specific energy outlays (in 
Joules/kg) on the range of the bound D is given in figure 9. 

In the supported phase of the bound from place, let the center 
of mass move with constant acceleration, while the distance 
traversed by the center of mass of the body, according to the 
coordinate c , is equal to Cnom"^min’ ^nom given 





nominal distance from the supporting surface to the center of 
mass of the body at the moment of lift-off from the supporting 
surface, and?f„iji is the given minimal permissible distance 
from the supporting surface to the center of mass of the body. 
For the examined vehicle, we will set Cnom=l»7 ni and ?mln=0.5 m, 
and we will determine the duration of the supported phase '^on* 
Then, the average specif ic forces per unit of weight, developed 
in the supported phase of the bound from place, are equal to 
v=i^Ton • The dependence of the average specific forces v ( in 
watts/kg) on the range of the bound is given in figure 10. 

We will study the question of the dependence of the 
specific energy outlays and the average specific outputs from 
acceleration of the force of gravity. Let 9, as before, be 
the acceleration of the force of gravity on Earth, and g* be 
the acceleration of the force of gravity on a planet on which 
the bounding vehicle is located. In view of t?ie fact that, 
with a fixed range of the bound the square of the velocity of 
the center of mass of the vehicle at the moment of lift-off from 
the supporting surface is proportional to the acceleration of 
the force of gravity, and the duration of the supported phase 
of motion is inversely proportional to the initial velocity of 
the center of mass, we obtain 



Thus, for example, for a vehicle moving along the surface of 
the moon, g*=l,62 m/sec. Consequently, on the ^surface of the 
moon, the specific energy outlays are roughly 6 times, and 
the average specific outputs are roughly 15 times, less than 
on Earth. The obtained results make it possible to draw a 
conclusion on the possibility of the effective utilization of 
bounding vehicles under the conditions of planets with a low 
force of gravity. 

The results obtained in the given paragraph can be /'l8 

generalized, utilizing the methods of similarity and dimen- 
sionality in mechanics [2]. With a proportional change in 
all of the linear dimensions by Z times, and the mass of the 
body and the leg sections by m times, the basic characteristics 
of motion change in the following manner: 

duration of unsupported and supported phases of the 
bound by VJI times, 

velocity of the center of mass of the vehicle and the 
body by VT times, 

angular velocity of the body and velocities at the leg 
joints by l//ji times. 


40 


energy outlays in the supported phase of the hound from 
place hy times, 

specific energy outlays V hy f, times, 
average specific outputs v hy times. 


Conclusion ^ 

Examined in the study is the problem of constructing 
the programmed motion of a hounding vehicle in the flight 
phase. According to the results of the study, one can draw 
the following conclusions: 

1. A mathematical model of the spatial motion of the 
hounding vehicle in the flight phase is constructed. ^The 
primary integrals of the equations of motion are obtained 
(law of motion of the center of mass of the vehicle and law 
of conservation of the moment of linear momentum relative to 
the center of mass) . 

2. A method is proposed for constructing- the translational 
motion of the legs in a relative system of coordinates, which 
ensures non- impact of lift-off and the required softness of 
placement of the legs on the supporting surface . 

3. An algorithm is constructed for constructing the_ 
programmed motion of the vehicle in the flight phase, which 
calculates the value of the linear and angular velocity of 
the body at the moment of lift-off from the supporting sur- 
face, necessary for realization of the impending flight phase. 

4. The construction of the programmed motion in the_ 
flight phase is carried out according to a complete or simpli- 
fied scheme of calculations, if this permits symmetry of 
motion of the vehicle. In the case of working of the algo- 
rithm according to a simplified scheme, a planar motion is 
actually examined instead of spatial motion of the vehicle. 

5 . The basic volume of the calculations for construction 
of the programmed motion is associated with the solution of the 
boundary value problem for determining the initial angular 
velocity of the body. Here, one must numerically integrate 
the differential equations. 

6. For standard motion conditions, which encompass a 
broad class of the most frequently encountered motion con- 
ditions of the vehicle, a table is compiled of the dependence 
of the initial angular velocity of the body on the parameters 

of the impending flight phase. The value of the initial angular 
velocity of the body is determined from this table. The volume 
of stored information is comparatively small. The value of 


the initial velocity of the center of mass of the body is 
computed according to finite formulas. Such an approach 
sharply reduces the time of calculation, and makes it possible 
to solve the problem of constructing the programmed motion for 
the flight phase in the process of motion of the vehicle. 

7 . The system of construction of the programmed motion of 
the bounding vehicle in the flight phase, which is part of its 
motion control system, automatically selects a complete or 
simplified scheme of calculations both for non-standard and 
for standard motion conditions. 

8. The algorithms for the system of construction of the 
programmed motion of the vehicle were worked out for the un- 
supported phase of motion by means of mathematical simulation 
on a computer. 
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Coefficients of the Integral of the Kinetic Moment 
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where 
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